THE SEMICLASSICAL EXPANSION OF THE T-J MODEL 
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The discovery of high temperature superconductivity has spurred intense investiga- 
tions of the two dimensional doped antiferromagnet. In the strong coupling limit, 
the t-J Hamiltonian, derived from the large-U Hubbard model, is often used to 
describe the low lying excitations. At zero doping, it directly reduces to the quan- 
tum antiferromagnetic Heisenberg model (QHM). Substantial progress has been 
recently achieved in understanding the Heisenberg limit, both theoretically and ex- 
perimentally [CI]. The effects of doping, however, are still highly controversial. As 
demonstrated in the other lectures of this Winter School, the problem of interacting 
spins and charges requires novel theoretical approaches. 

Let us start from the simplest one-band form of the Hubbard Model: 
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where c\ s creates an electron at site i with spin index s, and denotes a 

summation over all sites and their nearest neighbors on the square lattice. The 
large-U limit of Eq. (1.1) can be derived using second order perturbation theory 
in t/U. In the restricted Hilbert space of no double occupancies, the lowest order 
effective Hamiltonian is the t-J model: 
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where J = 4t 2 /U is the " super -exchange" constant, and P s projects onto the states 
with no double occupancies. (i]jk) are triads of nearest neighbors. At half fill- 
ing, (one electron per site), Eq. (2) reduces to the antiferromagnetic Quantum 
Heisenberg Model (QHM) on the square lattice: 
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where 

Sf =^E C L<'4' ' «= 1,2,3 (1.4) 



ss' 



are spin 1/2 operators in the subspace of one electron per site, and o a are the three 
Pauli matrices. Charge fluctuations are completely suppressed at half filling, and 
we are left with the spin excitations of the Heisenberg model. The two dimensional 
quantum Heisenberg model has been successfully addressed by two complemen- 
tary methods: the semiclassical approximation (large-S expansion) [CI], and the 
Schwinger boson mean field theory (large-N expansion) [Al]. There are also Fermion 
large-N expansions which do not recover the known properties of the QHM on the 
square lattice, and we shall not review them here. In order to address the doped 
system in a formulation than readily recovers the Heisenberg limit, we introduce 
two commuting Schwinger bosons aj,6j, and a slave fermion fa, to represent the 
allowed states of the projected Hilbert Space. The operators obey 



[a h a\] = , {fijj} = Sij , 
=0 , etc. 
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and satisfy the local constraints 

ajai + blbi + fj/i =1 Vi. (1.6) 
The t-J model is faithfully represented by 

H tJ = tJ2 flfjfij -j E - flfi) AtjA k j (1 -f}fj) , (1.7) 
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where 



4 = a \b)-b\a) 
T\j =a)ai + b)bi 



(1.8) 



It is now possible to generalize the t-J model to 5" > 1/2 by replacing the 
constraint (1.6) by 

alai + bjbi + fjfi =2S Vi. (1.9) 

In the half-filled (undoped) case, f'f =0, the Schwinger bosons describe states 
of spin S, and their bilinear forms yield the matrix elements of the standard spin 
operators 

^UV(^) =S? , a= 1,2,3. (1.10) 

In the presence of a hole at site i the spin size at that site gets reduced by 1/2. 
We shall see that the large-S limit yields a non trivial classical (S — > oo) theory 
provided we hold the scaled coupling constants J, i fixed, where 



J =AJS 2 , t =2tS 



(1.11) 



It is also possible to generalize the t-J model to large-N by introducing N flavors 
of bosons per site [A2]. The large-N mean field theory predicts spiral magnetic 
phases at finite doping concentrations [Jl]. Recently however, the fluctuation de- 
terminant has been found to be unstable (negative) in a range of momenta. The 
offending fluctuations were identified as local enhancements of the spiral distortion. 
Clearly, the holes drive strong perturbations of the spins on the lattice constant 
scale. These are difficult to treat by direct application of continuum and mean field 
approximations on the Hubbard and t-J models. This is why we shall concentrate 
on the large-S expansion. 

2. SPIN-HOLE COHERENT STATES 

Spin coherent states have been fruitfully used by Haldane to map the QHM onto 
the non linear sigma model [HI]. The construction of the spin coherent states path 
integral allowed a simple derivation of the topological Berry phases. The large-S 
expansion of this path integral provides a unified semiclassical treatment of the 
ordered and disordered phases of the quantum antiferromagnet. In the ordered 
phase, one obtains the usual spin wave expansion about the Neel state, which 
yields a very good approximation for the ground state staggered magnetization and 
uniform susceptibility even for spin 1/2 [CI]. 

In this lecture we generalize this approach to represent the t-J model by 
defining "spin-hole coherent states" . This will allow us to treat the short range in- 
teractions carefully, while observing the local constraints. We derive a semiclassical 
expansion of the ground state and low excitations in the presence of holes. Although 
the expansion is formally controlled by the large spin size S, at low doping we may 
rely on the success of this approximation for the S=l/2 QHM. 

Spin coherent states of spin S are defined as 

\Cl s (0, </>)> = (2S\)- 1 ' 2 (ua^ + vtf) 2S \0, 0) , (2.1) 
where O is a unit vector, and 

u = cos(6>/2)e^ /2 , v = sin(6>/2) e -^ /2 . (2.2) 

9, 4> are the lattitude and longitude angles on the unit sphere. 
The spin-hole coherent states are defined as follows: 

|n,0s = Ps) |o> + |n s -*> £/ f |o> • (2.3) 

where £ is an anticommuting Grassman variable. The states (2.3) allow a resolution 
of the identity, while conserving the constraint (1.9) for any given S-sector: 
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The factor as = (2S + 1)/(2S) is required for normalizing the matrix elements 
to unity. We shall be able to replace a by unity, by renormalizing the chemical 
potential [i in the grand cannonical ensemble. 



By repeatedly inserting the resolution of the identity (2.4) in the Trotter 
product expansion of the density matrix, we following the standard procedure [HI] 
and construct the path integral for the partition function: 

Z = JvCi vc Vi 



exp 
where 
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HtJ = (Am^ (2 . 6) 



is a function of the variables Oi,(f>i,£i, and 

A -S\ = (1- cos 9i(r)) far) (2.7) 

is the spin-kinetic term. The coordinate invariant notation uses the gauge potential 
A(O) which describes the vector potential of a unit magnetic monopole at the center 
of the sphere. The integral of the (2.7) over a closed loop on the sphere yields the 
solid angle enclosed by that loop. 

The Grassman "time derivative" is defined in its discrete form: 

e=(e(r)-e(r-e))/6, (2.8) 

where e is the infinitesimal timestep. 

H tJ has quadratic and quartic terms in the Grassmann variables. It is possible 
to integrate out the Grassmans exactly only for quadratic hamiltonians. Therefore, 
we decouple the four-fermion terms by a Hartree-Fock approximation. For our 
purposes this approximation will be justified by the following arguments: (i) hole- 
correlation corrections to this Hartree-Fock decoupling are higher order in hole 
density and (ii) the quartic terms vanish in the ferromagnetically correlated regions, 
where the hole density will turn out to be significant. 

Thus we arrive at a model which is quadratic in the Grassman variables: 
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7ij and iftij are the phases of u*Uj+v*Vj and (u*vl—v*ul)(ujVk—VjUk) respectively. 
The hole density matrix Pij[Cl] = (fjfj) is to be determined self-consistently, by 
solving for the Fermion ground state in the presence of a spin configuration {O}. 



The hole hamiltonian describes two distinct hopping processes: inter- 
sublattice hopping ("t-terms") and intra-sublattice hopping ("J-terms"). When 
the spins fij have short range antiferromagnetic order the 7 phases in the t-terms 
fluctuate wildly, while ipij = r]iA N ■ (xj — x^) represents a slowly varying Neel gauge 
field A N (x) which obeys 
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n ps r/itii is the staggered magnetization (Neel) field, where rji = +1, ( — 1) on 
sublattice A (B) is the "sublattice charge". V x is the "topological charge 
density" whose integral for continuos fields on a compact surface, is invariant [PI]. 

Weigmann, Wen, Shankar and Lee have studied Lagrangians which contain 
similar intra-sublattice A^-coupled hopping terms in the context of high T c super- 
conductivity [Wl]. Their starting model differed however from the t-J model in that 
it did not contain the inter-sublattice hopping terms (t-terms). These terms are 
not A N -gauge invariant, and do not conserve the sublattice charges. The t-terms 
cannot be neglected, especially in the t/J>l regime. We shall see, however, that 
in the large-S limit they are dynamically eliminated from the low energy and long 
wavelength Lagrangian. 

We begin by integrating out the fermions to obtain a purely spin partition 
function 
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(2.12) 



is the fermion free energy in the presence of a general spin history O(r). T T is 
the time ordering operator. Here we concentrate on the zero temperature case 
(3 = 00. Eq. (2.11) is a useful starting point for the semiclassical approximation. 
In the classical limit, S — > 00, the kinetic term is so large, that the only important 

configurations in the path integral are classical, i.e. the spins are frozen and (O) = 0. 
One has therefore to minimize H J + E* with respect to O for a given chemical 
potential or a specified number of holes. The second step is include the semiclassical 
fluctuations whose dynamics are controlled by the kinetic term. We discuss the 
single hole and the many hole cases seperately. 



3. RESULTS 



In Ref. [A3] , a Lanzcos algorithm on the Connection Machine was used to minimize 
the energy for 128 x 128 spins. The following results were found: 

The single hole: 

In the regime i/J> 0.87, we have found that the "polaron" which is a lo- 
cal alignement of spins, yields a lower energy than any of the possible uniform 
states, including: the Neel state, spiral states and canted states. This result helps 
to explain why the uniform spiral states were found to be unstable against short 
wavelength distortions in the mean field theory [A2]. The polaron variational pa- 
rameters were chosen to describe a ferromagnetic core, a transition region, and a far 
field antiferromagnetic tail. The latter is completely determined by the boundary 
condition 89 and the pure Heisenberg model (i.e. the Laplace equation). 

Our results are quite simple. For 1 < t/J < 4.1 the single hole energy is 
minimized by the five-site polaron (one flipped spin), depicted in Fig. 1. 



Fig. 1: The five site polaron. The hole density pu 
is primarily concentrated on the sites of the unfilled arrows. 
The circular arrows represent an allowed tunneling path, where 
the polaron hops two lattice constants to the left. 



The hole density is approximately 1/2 and 1 /8 on the central and neigboring 
sites respectively, with a small amount of leakage (due the J-terms in (2.9)) to sites 
further away. For 4.1 < t/J < 6.6, the polaron has two flipped spins (diagonally 
across a plaquette), and at larger values the core radius increases slowly as R c ~ 
(t/J) 1 / 4 , and the energy goes asymptotically as + At ~ (JT) 1 / 2 . The most 
important fact is that the small polarons do not distort the Neel background. In 
particular, the configurations centered on a bond [A4] are considerably higher in 
energy. We also find that the polarons have no tails [SI], i.e. 89 = 0, throughout 



the regime discussed above. This follows from competing contributions of order 
±J(59) 2 of H J and E* . Since in addition, the density p is exponentially localized 
near the polaron sites, we conclude that the classical interactions between polarons 
are short ranged. 

The polaron breaks the lattice translational symmetry. This symmetry is 
restored by tunneling events, where two spins i and k simultaneusly flip their ori- 
entation (see Fig. 1). The tunneling matrix element (the polaron's hopping 
rate), is non-perturbative in S~ x : 

T lk = r exp ^-^y* j « S 1 ' 2 ^ eMSa ik ) , (3.1) 

Eq. (3.1) can be computed as follows: The azimuthal coordinates are analytically 
continued — > <j>i = (j)' + i(p" , while their cannonical momenta Sf = (2S — pi) cos 9{ 
are kept real. It can be readily verified that ^^Sf, and H J + E? are conserved 

along the tunneling path S£(<f)) which minimizes the action. As a result of these 
conservation laws, we obtain a selection rule: tunneling can only take place between 
sites on the same sublattice! This, in effect, amounts to the elimination of the 
inter-sublattice "t-terms". 

otikjfiik are slowly varying dimensionless functions of t and J. For five site 
polarons, the number of spins involved in the tunneling path is at least three. For 
S = ^ we estimate the exponent to be roughly unity, but a fuller treatment of the 
multidimensional tunneling problem is required for a quantitative determination of 
r and the polaron's effective mass. 

The single polaron in a perfect Neel background occupies a Bloch wave of 
dispersion ek = 2T c (cos(2k x ) + cos(2k y )) + 2Tb(cos(k x + k y ) + cos(k x — k y )), where 
c, b denote the site of the other flipped spin as labelled in Fig. 1. By energetic 
arguments, r& < r c . A Berry's phase calculation of the tunneling matrix elements 
for S=l/2, yields an overall positive sign for r c ,rb > 0. Thus the single polaron 
energy is minimized at k = (n/2,n/2). This result agrees with other studies of 
the single hole spectral function in the t-J model [S2]. For small deviations of the 
background spins from antiferromagnetic order the tunneling rate Fik is modulated 
by the overlap of the background and the perfectly antiferromagnetic configurations. 
This overlap is just exp[irjiA N (xi — x^)]. and rji are the aforementioned Neel 
gauge field and sublattice charge respectively. We notice that A N couples in a 
gauge invariant way to the polarons, and that the sublattice charges are conserved 
in the hopping. 

Interactions : 

The interactions between two polarons were computed in the regime if J = 
1 — 4. We define Uf- = \eij — 2e^, where is the relaxed energy of a two-hole 
polaron with flipped spins at sites % and j. is repulsive, and of order 0.6 J — 2. 6 J. 
The intersite interactions, for neigboring polarons at sites (a)-(d) (see Fig. 1.) are 
plotted in Fig. 2. 

We find both attractive and repulsive interactions, and it is interesting to 
note that for t/J < 1.8 there is a near neigbour attraction of antiferromagnetically 



Fig. 2: Classical interactions between polarons, in units of J. 
Lines (a)-(d) represent the second polaron positions as labelled in Fig. 1. 
The solid line represents the relative condensation energy per hole 
of thehole-rich phase (see text). 

correlated spins. We also consider the possibility of polaron condensation into hole- 
rich domains [II]. The condensation energy per hole is determined by minimizing 
it with respect to the spin configuration, and the density. The spins in the hole- 
rich domains align ferromagnetically, and the energy per hole is given by e/ m = 

— At + A\J BJtw. This results coincides with that of Emery, Kivelson and Lin [II], 
except that their quantum correction factor B=0.584 is here set to 1/2. In Fig. 
2., the condensation energy Ae c = e/ m — is plotted. We find that it becomes 
negative at t/J = 2.7, above which phase-separation will occur for large S. 

Attractive interactions and negative classical condensation energies may yield 
charge density waves or superconductivity in the ground state of the quantum 
model. However, if realistic intersite Coulomb repulsions are added to the t-J model, 
the short range interactions may change sign. In particular, phase separation may 
be supressed, or pushed to higher values of t/J. 

The information given above allows us to write the effective Lagrangian for a 
dilute system of small polarons: 

+ Tike™*"-*** p\pk + X>5 PiPi PjPi 

{i;jk) ij 

(3.2) 

Eq. (3.2) is the most important result. C s p describes a two charge system of 
spinless Fermions pi with short range interactions Uf- coupled to Heisenberg spins. 
The formation of polarons can be viewed as a strong short- wavelength dressing of 
the original f-holes by the spins. As a consequence, the uncomfortable t-terms have 
been conveniently eliminated, and the effects of holes on the spin background is 
short ranged. This Lagrangian describes the low lying excitations of the so-called 
t'-J Hamiltonian 

n t>J =^J2S i -S j + Aj2 FikplpkAtAjk (3.3) 



A major advantage of the t'-J model over the t-J model is that in the small concen- 
tration limit 5 << 1, (3.2) is amenable to the continuum approximation. Following 
Haldane [HI] the spin interactions can be relaced by the (2+1) dimensional non lin- 
ear sigma model, with 5 dependent renormalized stiffness constant and spin wave 
velocity. The precise evaluation of the sigma model parameters for finite 8 is beyond 
the scope of this paper, but we expect that above some critical density 5 > S c the 
ground state is disordered [C2], i.e. a "spin liquid". In the massive "spin liquid" 
phase, Eq. (3.2) reduces to Wiegmann, Wen , Shankar and Lee's model [Wl]: 
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+ O(p^pp^p) . . . 



where m is the effective mass at k = (7r/2,7r/2), and the "electromagnetic" Neel- 
fields are F Vfl = (d^A^f — d u A^). k is the inverse spin correlation length, which is 
also the coupling constant of the gauge field. Previous analyses [Wl] have concluded 
that the ground state of (3.4) is most likely an RVB-type superconductor. Lee 
argued that the pairing is caused by two effects: (i) attraction between the opposite 
charges induced by the Neel gauge field, and (ii) suppression of coherent single 
particle propagation due to fluctuating Bohm-Aharonov phases, while the pairs 
(jA.pL) propagate as free bosons. Both (i) and (ii) are only valid in the magnetically 
disordered phase, a pleasing feature which agrees with the phase diagrams of the 
copper oxide superconductors. 

Aside from the mechanism of superconductivity, the small polaron theory 
could be checked numerically on finite lattice Monte-Carlo simulations, and experi- 
mentally in the copper-oxides and other doped antiferromagnets. For example: the 
polaron size can be estimated by NMR techniques [Ml] , and its internal excitations 
could be probed by optical absorption. In the frozen moments regime, one expects 
the polarons to exhibit conductivity typical of weakly localized semiconductors [C3] . 
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